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Abstract.  A non-collisional  kinetic  model  for  a non-neutral  plasma  in  a Perming  trap  is 
presented.  This  model  describes  the  evolution  of  the  z-integrated  distribution  function  of  the 
particles,  taking  into  account  the  three-dimensionality  of  the  problem.  The  general  properties  of 
the  model,  in  particular  the  conservation  laws,  are  studied.  The  model  is  also  related  to  the  fluid 
model  proposed  by  Finn  et  al.  and  refined  by  Coppa  et  al.  Finally,  numerical  investigations  are 
presented  concerning  the  equilibrium  solutions  of  the  model. 

INTRODUCTION 

Usually,  the  dynamics  of  non-neutral  plasmas  confined  in  Penning-Malmberg  traps 
is  described  by  employing  a two-dimensional  drift-Poisson  model,  where  charged 
particles  are  regarded  as  straight  lines  (strings)  of  uniform  density,  due  to  the  very 
high  value  of  the  axial  bouncing  frequency.  According  to  this  model,  the  particles 
exhibit  an  ExB  drift  motion  and  the  electric  field  is  computed  self-consistently  from 
the  charge  density  by  using  the  Poisson  equation. 

This  model  predicts  possible  instabilities  for  the  azimuthal  modes  ingPl.  For  the 
mg=\  mode,  discrete  modes  are  always  stable  for  any  density  profiles  [1],  while  the 
continuum  spectrum  can  only  produce  an  algebraic  growth  [2].  Nevertheless, 
experiments  show  that  the  mode  growth  is  actually  exponential;  the  contradiction 
between  experimental  results  and  the  linear  two-dimensional  theory  is  a challenging 
problem  in  the  theory  of  non-neutral  plasmas. 

The  present  work  grounds  on  recent  studies  on  the  evolution  of  non-neutral  plasmas 
in  a Penning-Malmberg  trap,  pointing  out  the  important  role  of  kinetic  effects  and  of 
the  finite  length  of  the  device,  in  particular  for  the  ino=\  diocotron  instability  [3-5].  In 
fact,  when  a particle  approaches  the  border  of  the  plasma,  it  feels  a confining  potential 
which,  in  general,  depends  upon  both  the  radial  and  the  axial  coordinates.  The  radial 
component  of  the  confining  electric  field  causes  an  ExB  drift  in  the  azimuthal 
direction,  affecting  the  rotation  frequency.  The  drift  depends  upon  the  axial  energy  of 
the  particles,  which  affects  the  penetration  in  the  confining  potential  [6]. 

The  aim  of  the  present  work  is  to  develop  a kinetic  theory  for  non-neutral  plasmas 
taking  into  account  self-consistently  all  these  effects.  The  kinetic  model  assumes  that 
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the  frequency  are  ordered  as  fir » fi/, » fif,  being  fir  the  cyclotron  frequency,  fi* 
the  bouncing  frequency  and  fif  the  ExB  drift  rotation,  which  is  comparable  to  the 
frequency  of  the  modes  of  interest.  Within  this  assumption,  the  planar  motion  of  the 
electrons  is  well  described  by  the  dynamics  of  their  guiding  centers,  given  by  the  ExB 
drift  velocity,  and  the  complete  kinetic  description  of  the  plasma  is  then  provided  by 
the  distribution  function  f{r,9,  z,  t),  where  ^is  the  axial  energy  of  the  particles.  The 
description  is  completed  by  the  self-consistent  Poisson  equation  for  the  electrostatic 
potential  z,  /). 

In  order  to  reduce  the  complexity  of  the  problem,  the  dimensionality  of  the  model 
is  reduced  by  integrating  along  the  axial  direction.  If  the  electrostatic  potential  “seen” 
by  a string  of  electrons  varies  slowly  in  time,  the  ergodic  distribution  in  the  phase- 
space  (z,  Vj)  can  be  assumed  for  the  particles  of  the  same  energy  and  a kinetic  equation 
for  the  z-integrated  electron  distribution,  F{r,  9,  t),  can  be  deduced.  The  properties 

of  such  kinetic  equation  are  the  main  subject  of  the  present  work. 


DEDUCTION  OF  THE  KINETIC  EQUATION 

The  starting  point  in  deducting  the  present  model  is  the  kinetic  equation  for  the 
complete  distribution  function,/ {r,9,  z,  t). 


8t  " 
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8z^ ’ d4\dtj 


= 0 


(1) 


where  the  ExB  drift  approximation  is  used  and  the  phase-space  variable  ^ is  the  axial 
energy  of  the  particles 


^='^mv]-e<l>{e,e,z,t)  (2) 

while  the  time  derivative  of  ^is  simply  given  by  d^l  dt  = - ed(j>l  dt. 

The  z-integrated  distribution  F{r,9,^,t)  is  introduced  by  factorizing  the  distribution 
function  as 


f{r,9,^,z,t)  = F(r,9,4,t)g{r,9,^,z,t)  (3) 

The  function  g is  defined  assuming  the  particles  of  a given  energy  to  be 
distributed,  in  the  phase-space  (z,  v,),  according  to  the  ergodic  distribution.  Moreover, 
g is  normalized  in  order  that  f g(r,9,^,  z,  t)  dz  =1 . This  leads  to 


g(r,9,4,z,t) 


[e^{r,9,z,t)-^Y'- 
\[e4>(r,9,z,t)-Y\'''  dz 


(4) 
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Within  this  assumption,  the  following  equation  for  the  function  F{r,0,  t) 


+ (5) 

is  obtained  by  integrating  Eq.  (1)  along  the  axial  direction  and  defining  the  bounce- 
averaged  streaming  velocities  as 

(■e,  td<t> 

Vfl  = I — - — gdz  (6) 

The  kinetic  model  is  then  completed  by  the  self-consistent  Poisson  equation  in  the 
Pemring  trap 


VV  = — [ F(r,0,4,t)g(r,0,z,4,t)d^  (7) 

Finally,  it  must  be  observed  that  the  kinetic  equation,  obtained  using  the  adiabatic 
invariant 


m = + (8) 

as  phase-space  variable  instead  of  is  simpler,  in  principle;  however,  the  inversion  of 
Eq.  (8),  required  in  the  solution  of  the  Poisson  equation,  is  a challenging  numerical 
problem. 


GENERAL  PROPERTIES 

The  distribution  function  F can  be  related  to  physical  quantities  that  can  be 
measured  experimentally.  Denoting  with  o(r,0,  l)  the  z-integrated  plasma  density  (an 
experimentally-measurable  quantity)  and  with  F^(r,0,  t)  the  normalized  distribution 
of  the  kinetic  energy  of  the  particles,  4)  the  distribution  function  F can  be  expressed  as 
F{r, 0,^,1)  =o{r,0,  t)F^{r,0,  t),  with  being  (f>,„  the  potential  in  the  point 

where  the  kinetic  energy  of  the  particles  is  measured  and  | F^{(^,  r,0,  t)d<^=\. 

The  kinetic  model  presented  here  satisfies  some  conservation  laws.  The  canonical 
angular  momentum,  L defined  as 

- Y n 1 r^nr,0,4,t)drd0d^  (9) 

and  the  angular  momentum,  4,  proportional  to  the  axial  angular  momentum  of  an 
incompressible  fluid. 
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/,  = J e.  - r x(Ex  B)Jr 


(10) 


are  both  conserved. 

As  far  as  the  energy  is  concerned,  the  total  energy  of  the  plasma,  E,oi,  can  be  written 
as  the  sum  of  potential  and  kinetic  energies,  E ,„i=  E poi+  E kin,  being 


Ep„,  \e(pF{r,e,4,t)g(r,9,z,^,t)clrd4 

Eki„  = ^J  \mv:{r,e,4,t)g{r,9,z,^,t)drd4 


(11) 


The  calculation  of  r/E,oi  / dt  leads  to  the  conclusion  that  the  model  conserves  the 
total  energy,  as  the  variation  of  the  energy  of  the  plasma  equals  the  power  provided  by 
the  electrodes: 


= f^V^-nr/5  (12) 

dt  2 dt  i 

CONNECTION  WITH  A RECENTLY-PROPOSED  FLUID 

MODEL 

The  present  kinetic  model  can  be  related  to  the  fluid  model  proposed  by  Finn  et  al. 
[5]  and  refined  by  Coppa  et  al.  [4].  The  main  assumption  of  that  model  is  that  the 
Maxwell-Boitzmann  distribution  is  reached  along  the  axial  direction,  so  that  the 
particle  density  n(r,9,  z,  t)  can  be  expressed  as 


n(r,9,z,t)  = N{r,9,t)exp 


e(j>{r,9,z,t) 


kj 


(13) 


where  N{r,9,  t)  is  a function  that  does  not  depend  upon  z.  Within  this  assumption,  the 
continuity  equation  for  the  z-integrated  density,  a{r,9,  t),  is  written  as 


^ + V,(V,a)  = 0 


(14) 


where 


f V jex^[{etf)l{kj)\lz 

J exp[(e^J)/(^a7’)]r/z 


(15) 
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The  fluid  model  can  be  deduced  from  the  present  kinetic  model  by  integrating  Eq. 
(5)  with  respect  to  assuming  that  the  function  can  be  written  as: 


F{r,0,^,t)  = Fir,d,t)exp\ 


k,T 


dz 


yle(f>  + ^ 


(16) 


This  assumption  is  equivalent  to  affirm  that  the  electrons  are  distributed  according 
to  the  canonical  distribution;  i.e.  the  electron  density  is  proportional  to 


NUMERICAL  RESULTS 


The  electrostatic  potentials,  corresponding  to  axially-symmetric 

distributions,  F{r,^,  have  been  calculated  by  solving  numerically  the  self-consistent 
Poisson  equation  [Eq.  (7)].  The  equation  has  been  discretized  in  space  (using  a 
uniform  grid)  and  energy.  Then,  suitable  quadrature  formulas  have  been  used  to 
evaluate  the  ^ and  z-integrals.  The  Poisson  equation  has  been  solved  by  adding  a term 
B^dt  to  the  left  hand-side  and  seeking  the  steady  state  solution  of  the  diffusion-like 
equation  so  obtained. 
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FIGURE  1.  Plasma  density  difference,  with  respect  to  the  Maxwellian  distribution,  for  the  Maxwellian 
distribution  truncated  at  the  energy  I=ksTH  (dotted),  7=^^? (dashed)  and  l=2kBT  (solid).  The  confining 
potential  of  the  trap  is  F=-50  V,  the  central  electrode  length  Lc=26  cm,  the  end  electrodes  length  1^=4 
cm,  the  wall  radius  7?„,=3.5  cm,  the  plasma  temperature  r=1000  K,  the  magnetic  field  Bo=l  T and,  the 
parameters  of  the  density  profiles  introduced  by  Finn  [5]  are:  «o=5T0'’  m'^,  //=5,  r^=2  cm. 
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FIGURE  2.  Bounce-averaged  drifi  frequency  for  electrons  with  Maxwellian  distribution  at  different 
axial  energy:  ^ =20.5  eV  (solid),  ^ -23.3  eV  (dotted)  and  4 “26  eV  (dashed)  for  two  different  values  of 
L^.  The  predictions  of  the  2-D  theory  are  shown  in  thick  .solid  lino.  The  other  trap  parameters  are  the 
same  as  Fig.  1. 

The  equilibrium  solutions  have  been  studied  having  fixed  the  plasma  density  at  the 
center  of  the  trap  as  the  density  profiles  introduced  by  Finn  et  cil.  [5],  The  velocity 
distributions  have  been  chosen  to  be  truncated  Maxwellians,  as  suggested  by 
Hilsabeck  et  at.  [3]. 

Some  numerical  results  are  presented  on  Figs.  1-2.  Figure  1 points  out  the  effect  of 
the  velocity  distribution  on  the  plasma  density  by  comparing  different  truncated 
Maxwellians.  Figure  2 shows  how  the  rotation  drift  frequency  depends  on  the  particle 
energy. 
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